Deterministically entangling distant nitrogen-vacancy centers by a 

nanomechanical cantilever 
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We present a practical scheme by global addressing to deterministically entangle negatively 
charged nitrogen-vacancy (N-V) centers in distant diamonds using a nano-mechanical cantilever 
with the magnetic tips strongly coupled to the N-V electron spins. Symmetric Dicke states are gen- 
erated as an example, and the experimental feasibility and challenge of our scheme are discussed. 
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Over past years, negatively charged nitrogen-vacancy 
(N-V) centers in diamond have been considered as one of 
the most promising candidates for solid-state quantum 
information processing (QIP) due to their long coherence 
time at room temperature [l| . The intensive experimen- 
tal studies have been paid on quantum logical gating Q , 
state storage and transfer [3], and preparation of entan- 
gled states 4J, using the electron-spin and the nearby 
nuclear spins within a single diamond. 

In the context of a large-scale diamond-based QIP, the 
generation of multipartite entangled states among dis- 
tant N-V centers is of particular importance and exper- 
imental challenge. One of the proposed approaches is 
to use the measurement-based QIP scheme By de- 
tecting the emitted photons which have been entangled 
with different N-V electron spins, distant N-V centers 
can be entangled in a nondcterministic way. Such an 
approach has the major advantage that no direct inter- 
actions among qubits are required, which implies a rela- 
tively easy control. An alternative is the scalable optical 
coupling between N-V centers and high-Q microcavities 
. By sharing the common modes of a microcavity or 
sharing flying single-photon pulses, the remote N-V cen- 
ters could be coherently manipulated or entangled, which 
may open up the possibilities of N-V-based scalable QIP. 

To have an efficient processing of quantum informa- 
tion, however, we prefer deterministic manipulations. In 
the present work, we will show the possibility to entangle 
multiple N-V centers by a nano-mechanical cantilever in 
a deterministic fashion. We had noticed a recent work 
for a strong magnetic coupling between electron-spin de- 
grees of freedom of a N-V center and a nano-mechanical 
resonator 0, [1] • In our case, we have a line of N-V cen- 
ters, each of which is located in a nanoscale diamond. We 
will show how to generate the symmetric Dicke states Q , 
denoted by \D^'^) with equally weighted superposition 
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FIG. 1: (Color online) (a) Schematic of the physical system 
consisting of an array of N-V centers strongly coupled to the 
quantized motion of a cantilever with length I, width w and 
thickness r. An external magnetic field ^ext is applied along 
the symmetry axis of the N-V centers, (b) The energy level 
diagram for the ground triplet state of the electron spin within 
a N-V center, with a zero-field splitting Dgs between spin 
levels ms = and rus = ±1. We restrict our study to the 
down state ms = and the up state nis = — 1 by driving a 
microwave field. 



of all permutations of m excitations among n particles, 
by global addressing. Such multipartite entangled states 
own great properties of robustness against decoherence 
and particle loss and have been observed experimentally 
with linear optics systems • The favorable features of 
our scheme include the deterministic preparation of en- 
tangled electron-spin states between N-V centers and the 
potential to be scaled up with near future technologies. 
Besides, as a scheme with global addressing, only a mi- 
crowave field is required during the entanglement prepa- 
ration, which greatly reduces the experimental time and 
difficulty. Furthermore, our proposed scheme could be 
realized with very high fidelity due to our fast manipula- 
tions and the long decoherence time of electron spins in 
N-V centers. 
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As sketched in Fig. 1(a), an array of evenly spaced 
magnetic tips are attached at the end of a nano- 
mechanical cantilever, under which are n one-to-one cor- 
respondent N-V centers. In order to entangle the electron 
spins in spatially separate N-V centers by a cantilever, 
we have to couple the electron spins to the motion of the 
resonator. This can be achieved by exposing the elec- 
tron spins to the magnetic field gradient created by the 
magnetic tips. For clarity of our description, we label 
\j = XeB^gZQ {j — 1,2,- --n) as the coupling strength, 

where — \dz'^tipj\, Zq is the zero-point amplitude of 
the cantilever regarding the jth tip and Xe is the charge- 
mass ratio of an electron. For simplicity, we assume that 
all the N-V centers are equally coupled to the cantilever 
(Aj — A), which could be achieved by elaborately ad- 
justing the magnetic field gradient of the tips and the 
distances between the tips and the corresponding N-V 
centers. Furthermore, in most part of the paper, we will 
neglect the influence from the nuclear spins, which im- 
plies a pure N-V center without ^^C doped or with well 
polarized nuclear spins of ^^C or ^^N [ll|. For our pur- 
pose, we encode qubits in two sublevels: nis — and 
771s = — 1 of the ground triplet state, and a microwave 
field will be employed to drive the two sublevels. As a 
result, the Hamiltonian of the whole system can be writ- 
ten in units of /i = 1 as 

H= E:Li^o\-l)j{-l\ + 'ya^a 
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where ujo « Dgs — Xe\^ext\ is the energy splitting be- 
tween sublevels |0) and | — 1) under the externally ap- 
plied magnetic field Bext along z axis, i.e., the quanti- 
zation axis, Dgs (around 2.87 GHz is the zero-field 
splitting between the spin levels 777s = and ttIs = ±1 
shown in Fig. 1(b). i' and a (a^) are the frequency and 
annihilation (creation) operator of the motional mode of 
the cantilever. ^j(t) and Uj are the Rabi frequency and 
the frequency of the microwave regarding the cantilever 
couphng to the jth N-V center, and ai — \ — l)j (— 1| — 
\0)j (0|. For simplicity, we may assume in following treat- 
ment ^j{t) = Q{t) and ujj = w (j = 1,2, • • -77). So the 
Hamiltonian of the total system in the interaction picture 
is 
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If w = Wo, Eq. (2) in a new basis |±) 
reads as J^^^^i + E"=i 

= 1+).,- (-1 + I-).- (+1 and ai = |+), (+| - h), (-| . 
Considering il{t) is comparable to ly, we rewrite Eq. (2) 
in a rotating frame under rotating wave approximation. 



obtaining 



Hr = 



En A 
.1=1 2 



|exp 



d{t)dt 



ti 



H.c. 



(3) 

with 5{t) = v- n{t). 

From the viewpoint of QIP, the vibrational mode of a 
cantilever could be considered as a common data bus for 
distant N-V centers. Therefore, using Eq. (3), we may 
achieve various kinds of multiple entangled states. As an 
example, we will show below how to prepare symmetric 
Dicke states of the electron spins in distant N-V centers 
by global microwave addressing. 

We consider a Si cantilever beam of length I — 5 /im, 
width w — 50 nm and thickness r = 50 nm with effective 
mass Af ~ 7.28 x 10^^* kg, resonant frequency of the 
lowest-order flcxural mode ~ 14.87 MHz and the zero- 
point amphtude of fluctuation ^l/(2Mi/) ~ 6.98 x lO^^'^ 
m 13]. The N-V-cantilever coupling strength A ~ 0.96 
MHz has been achieved under a magnetic field gradient 
7.8 X 10*^ Tm-i at a distance ~ 25 nm away from the 
tip 0, addition, the distance between the nearest- 

neighbor tips is set to be around 150 nm. 

We will generate two nontrivial three- 
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after performing a time-dependent phase transformation. 
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FIG. 2: (Color online) Left: The energy level pattern of a 

string of three NV centers with an initial state | 

globally driven by a microwave field, where we set 5{t) = 
and Dp is for the motional state of the cantilever. Right: 
The MS transformation produces a two-level subsystem where 

|/3i) = \D[^^), \a) = \ ) and two decoupled dark states 

Ift) = (I + ) - 2| - +-) + I - -+))/V6 and = (| + 

— )-\--+))/V2. 



Let us first consider how to prepare \D 



We 



sume that the vibrational state of the cantilever has been 
cooled to the Fock state |l)p 0, [sl and the electron-spin 
states of N-V centers are initially prepared in 



By means of Morris-Shore (MS) transformation [15|, it is 
possible to factorize the Hamiltonian into a smaller closed 
subspace [l^. As shown in Fig. 2, in MS basis the orig- 
inal system is decomposed to a two-level subsystem plus 
two independent single-level subsystems. The effective 
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Hamiltonian of the two-level subsystem is written as 
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(4) 



in the basis {|/3i)|0)p, |a)|l)p}, where = and 

|a) = I ). The top state |/3i) is connected to the 

initial state \a) with an effective coupling strength \/3X/2 
and can be interchanged by a microwave field addressing 
the three N-V centers globally. When 6{t) is changed to 
zero (by using a nearly square pulse of the microwave), 

the desired loi^'^) can be realized at t=7r/ {VSX) with unit 
probability. 
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FIG. 3: (Color online) (a) Left: The energy level pattern of 

a string of three NV centers with an initial state | )|2)p 

globally driven by a microwave field. Right: The MS transfor- 
mation produces a three- level subsystem where l^i) — [D^^^), 

|?7i) = \D^''), Ix) = I ) B^nd two decoupled two-level 

subsystems 1^2) = (| + +-) - 2| H h) + | - ++))/V6, 

\V2) = (I + — ) - 2| - +-) + I - -+))/V6 and = 
(I + +-> - I - ++))/V2, \m} = (I + — ) - I - -+})/V2. 
(b) Adiabatic evolution of the three-level subsystem, where 
the three black-dotted lines represent and |x) re- 

spectively. The up state and ground state \x} can be 
interchanged by adiabatic passages (arrows) under the adia- 
batic condition, (c) Populations of states |7;i), \x) during 
the adiabatic evolution. 
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Now we turn to the preparation of ').ln contrast to 

|I?3^'') preparation, the state of the whole system in this 

case is initially prepared in | )|2)p. After performing 

MS transformation, we obtain a three-level subsystem 
plus two additional two- level subsystems [See Fig. 3(a)]. 
The effective Hamiltonian reads as 
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in the basis {|Ci)|0)p, |?7i)|l)p, |x)|2)p} , where 



\D 
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D^3^) and Ix) 



It is easy to 



check by setting S{t) = that the population proba- 
bility of |I?3 ') will only reach 0.95. In order to gener- 

(2) 

ate |Z?3 ) with higher probability, we will employ adi- 
abatic techniques as below. The top state and 
the bottom state |x) in Fig. 3(b) can be interchanged 
through adiabatic passages under the adiabatic condi- 

< 1 



{Ei{t)\dtH^ff\m))/[Ek{t)-Ei{t)r' 



tion: max 

{I ^ k), where \Ek) and Ek are the instantaneous eigen- 
states and eigeN- Values of Eq. (5). In Fig. 3(c), the 
populations of jryi) and |x) are numerically calcu- 
lated from Eq. (5), where we have assumed a purely 
linear microwave chirp (using sawtooth pulse of the mi- 
crowave). By addressing the three N-V centers globally, 

('21 

we could achieve the population probability of jUg ) to 
be higher than 0.99 through this adiabatic passage. 

With similar steps to the above cases for three N-V 
centers, a generalization to \D^'^) with 7i > m > 3 is 
straightforward. 

In the above treatment, we have assumed the N-V cen- 
ters to be equally coupled to the cantilever. In realistic 
case, however, the coupling should be slightly different 
due to imprecision. Figure 4(a) demonstrates a simpli- 
fied case for the fidelities of the produced states under 
the relative error of the coupling strength A, defined by 
A„ = [1 — (n— 1) A]A, where the fidelities still remain high 
in the case of the small relative error. 

The main decoherence sources associated with our im- 
plementation are dephasing of the electron spins and 
the intrinsic heating of the cantilever. Since our op- 
eration for the symmetric Dicke states could be ac- 
complished on the time scale of 1.8 ^ 30 /iS, which 
is much shorter than the reported electron-spin relax- 
ation time Ti = 6 ms Q as well as the dephasing time 
T2 — 350 fjLS induced by the nuclear-spin fluctua- 
tion inside a N-V center. Thus, the influence on the 
N-V electron-spin from the intrinsic damping and de- 
phasing is negligible. On the other hand, as the can- 
tilever is ultra-sensitive to the magnetic force generated 
between the tips and the N-V electron spins [l3|, [3, [3] , 
we will concentrate on the heating effect of the can- 
tilever by the master equation dtp = —i[H^^j,p]-\-hp 
with k =1, 2 [corresponding to Eqs. (4) and (5), respec- 
tively], where Lp = r(fi.t;i-|-l) {2apa^ — a^ap — pa^aj /2+ 



l/[exp(^) 



Tnth{2a) pa- aa^ p- paa^) /2, nth 

and ks is the Boltzmann constant. In general, the heat- 
ing factor of the cantilever is described hy T = v/Q [l3[, 
where u and Q are the vibrational frequency and the 
quality factor of a cantilever. Numerical calculations are 
performed for the fidelities of the Dicke states in Fig. 



4(b), where the fidelities of \D^^'') and \D'^') could reach 
0.998 and 0.996 respectively in the case of the Q value 
10^. If the Q value is lowered to 10^, the fidelity drops 



l(2) 



to 0.985 for jcl^h and to 0.957 for \D. 
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FIG. 4: (Color online) The fidelity of symmetric Dicke states l-Da^') (solid line) and \Di/^) (dashed line) versus (a) relative 
error A of the coupling strength; (b) heating factor F of the cantilever. 



We address some remarks for experimental issues of 
our scheme. To make sure that the N-V centers are all 
strongly and nearly equally coupled to the cantilever, 
we have to restrict the tips within a small region near 
the head of the cantilever, which implies a limitation of 
the number for available N-V centers (at most tens of 
N-V centers in our present design). So a potential for 
scalability with our scheme might be to combine with 
the non-deterministic method [5]. Besides, we require 
nanoscale diamonds with a single N-V center in each and 
the N-V centers should be located very close to the sur- 
face of the diamonds. To our knowledge, such diamond 
nanoparticles are not in principle unavailable although 
they have not yet been achieved experimentally. On the 
other hand, due to the small spacing between adjacent 
N-V centers, individual addressing with a microwave field 
is of difficulty. As a result, although entanglement could 
be achieved with global operations, to carry out universal 
quantum computing, we may have to employ the near- 
field optical technique [T^ with spatial resolution of the 
order of tens of nanometers, which had been widely used 
in optics and surface science, and recently also employed 
for QIP Alternatively, a magnetic field gradient 

along X-axis would help to distinguish different qubits in 
frequency space. Thus universal QIP is available, and 



robust quantum gating, such as composite pulse method 
[2H or bang-bang control [22] , employed previously could 
be straightforwardly applied to our scheme. Finally, cool- 
ing a cantilever with many magnetic tips to a required 
low temperature (114 fiK in our design) could be of great 
technical difficulty with current technology but rnay be 
resolved by recently presented cooling schemes 0, Q- 

Further consideration would iN-Volve nuclear spins 
which are good at storing quantum information due to 
the long decoherence time 23]. In this sense, we may 
also consider to encode the qubits in the nuclear spins, 
and employ the electron spins as ancillas. The hyperfine 
interaction helps to transfer the generated entanglement 
of electron spins to the corresponding nuclear spins of 
N-V centers [3]. 

To conclude, we have proposed a scheme to deter- 
ministically entangle a string of N-V centers by strongly 
coupling to the quantized motion of a nano-mechanical 
cantilever. Due to the long relaxation and dephasing 
time of the electron spins in N-V centers together with 
the weak intrinsic heating effects of the cantilever, our 
present scheme would enable a promising way for future 
N-V-center-based QIP. 

The work is supported by NNSFC under Grant No. 
10774042, by CAS, and by NFRPC. 
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